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Abstract 

Spatial distributions of energy deposited by an extensive air shower in the atmo- 
sphere through ionization, as obtained from the CORSIKA simulation program, are 
used to find the fluorescence light distribution in the optical image of the shower. 
The shower image derived in this way is somewhat smaller than that obtained from 
the NKG lateral distribution of particles in the shower. The size of the image shows 
a small dependence on the primary particle type. 
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1 Introduction 



The fluorescence method of extensive air shower detection is based on record- 
ing hght emitted by air molecules excited by charged particles of the shower. 
For very high energies of the primary particle, enough fluorescence light is 
produced by the large number of secondaries in the cascading process so that 
the shower can be recorded from a distance of many kilometers by an ap- 
propriate optical detector system [1,2]. As the amount of fluorescence hght is 
closely correlated to the particle content of a shower, it provides a calorimetric 
measure of the primary energy. 

The particles in an air shower are strongly coUimated around the shower axis. 
Most of them are spread at distances smaller than several tens of meters from 
the axis, so that when viewed from a large distance, the shower resembles 
a luminous point on the sky. Therefore, a one-dimensional approximation of 
the shower as being a point source might be adequate in many cases regard- 
ing the shower reconstruction. For more detailed studies, however, the spatial 
spread of particles in the shower has to be taken into account. This is espe- 
cially important for nearby showers, where the shower image, i.e. the angular 
distribution of hght recorded by a fluorescence detector (FD), may be larger 
than the detector resolution. 

The image of a shower has been studied in Ref. [3] , where it was shown that for 
a disk-like distribution of the light emitted around the shower axis, the shower 
image has a circular shape, even when viewed perpendicular to the shower axis. 
Analytical studies including lateral particle distributions parameterized by the 
Nishimura-Kamata-Greisen (NKG) function or estimates based on average 
particle distributions taken from CORSIKA [4] were discussed in Ref. [5] and 
Ref. [6], respectively. 
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In this paper, detailed Monte Carlo simulations of the shower image based on 
the spatial energy deposit distributions of individual showers are performed. 
By using the energy deposit of the shower particles as calculated by COR- 
SIKA [7], the previous simplified assumption of a constant fluorescence yield 
per particle is avoided. Assuming a proportionality between the fluorescence 
yield and ionization density, the light emitted by each segment of the shower is 
determined. A concept is developed to treat the shower as a three-dimensional 
object, additionally taking into account the time information on photons ar- 
riving at the FD. In contrast to previous analytical studies, shower fluctua- 
tions as predicted by the shower simulation code arc preserved and studied. 
Propagation of the light towards the detector, including light attenuation and 
scattering in the atmosphere is simulated, so that the photon flux at the detec- 
tor is calculated. The resulting distribution of photons arriving simultaneously 
at the detector, i.e. the shower image, is compared to results obtained by using 
the NKG approximation of particle distribution in the shower. The compari- 
son is performed for different shower energies and different primary particles. 
In particular, it is checked whether the shower width depends on the primary 
particle type. 

The plan of the paper is the following: definition of the shower width and 
algorithms of fluorescence light production are described in Section 2. In Sec- 
tion 3 the size of shower image in the NKG and CORSIKA approaches is 
calculated and its dependence on primary energy, zenith angle and primary 
particle is discussed. Conclusions are given in Section 4. 
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2 Simulations 



2. 1 Shower width and shape function 

Given an optical imaging system for recording the light emitted by a shower, 
the shower width is defined as the minimum angular diameter a of the image 
spot containing a certain fraction F{a) of the total fight recorded by the FD. 
The image is considered to be recorded instantaneously, i.e. with an integration 
time such that the corresponding angular shower movement is well below the 
angular resolution of the detector. 

Four main components of light contribution can be distinguished: (i) fiuores- 
cence light, with isotropic emission; (ii) direct Cherenkov radiation, emitted 
primarily in the forward direction; (iii) Rayleigh-scattered Cherenkov light; 
(iv) Mie-scattered Cherenkov light. The relative contributions of these compo- 
nents depend on the geometry of the shower with respect to the detector [8] , 
but in most cases the fluorescence fight dominates the recorded signal. As- 
suming only minor effects on the shower width by absorption and scattering 
processes during the fluorescence light propagation from the shower to the 
detector, the light fraction F{a) is mainly determined by the corresponding 
light fraction F{r) emitted around the shower axis 



where /(r) is the (normalized) lateral distribution of fluorescence light emit- 
ted. The main task is therefore to derive /(r), which is also referred to as the 
shape function, since the brightness distribution of the shower image depends 
on the shape of /(r). The shape functions in different methods of evaluating 
fluorescence fight production described in the following, i.e. in the NKG and 
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CORSIKA approach, and for different primary particle types in the CORSIKA 
approach will be compared. 

Photon propagation towards the detector is simulated based on the Hybrid_f adc 
simulation software [9], including Rayleigh scattering on air molecules and Mie 
scattering on aerosols. The final shower image is constructed by recording the 
photons that arrive simultaneously at the detector [5]. These photons that 
form an instantaneous image of the shower, originate from a range of shower 
development stages. Thus, for a precise description of the shower image, we 
need to take into account also the geometrical time delays of the photons 
coming from these stages, as will be discussed later. 

Since this work is intended as a general study, the resulting photon distribu- 
tion after light propagation is assumed to be recorded by an ideal detector. 
Possible effects of specific detector conditions such as spatial resolution or 
trigger thresholds will also be commented on, however. Investigations specific 
to the fluorescence detectors of the Pierre Auger Observatory are described in 
Ref. [10] 

2.2 Fluorescence light production 

As the shower develops in the atmosphere, it dissipates most of its energy 
by exciting and ionizing air molecules along its path. From de-excitation, UV 
radiation is emitted with a spectrum peaked between 300 and 400 nm (three 
major lines at 337.1 nm, 357.7 nm, 391.4 nm). Measurements have shown that 
the variation of the fluorescence yield n^^o, i-c the number of photons emitted 
per unit length along a charged particle track, as a function of altitude is quite 
small for electrons of constant energy. For example, the measured fluorescence 
yield of an 80 MeV electron varies by less than 12% around an average value of 
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4.8 photons/m over an effective altitude range of 20 km in the atmosphere [11]. 
This motivates to some extent the use of a constant, average fluorescence yield 
per shower particle, as will be described in the NKG approach (section 2.2.1). 

On the other hand, since the fluorescence light is induced by ionizing and 
exciting the molecules of the ambient air, the fluorescence yield is expected 
to depend on the ionization density along a charged particle track [11,12,13]. 
Most shower particles contributing to the energy deposit in air have kinetic 
energies from sub-MeV up to several hundred MeV [7] which is in the energy 
range of considerable dependence of ionization density on particle energy. As 
an example, a measurement of the air fluorescence yield [11] between 300 and 
400 nm at pressure 760 mm Hg is shown in Figure 1. The sohd hne repre- 
sents the electron energy loss dE/dX as a function of the electron energy. 
The minimum of this curve corresponds to 1.4 MeV electrons with energy 
loss {dE / dX)\iAMeV — 1-668 MeV/gcm~^ and fluorescence yield n^,o = 3.25 
photons per meter. We note that dE/dX increases by about 50% for energies 
from 1.0 MeV to 100 MeV, so the energy spectrum of electrons in a shower 
and its variations with atmospheric depth should be taken into account for 
an accurate determination of the fluorescence emission of the shower. There- 
fore to obtain a more realistic simulation of the spatial distribution of light 
production, the distribution of the energy deposit in the shower is used in the 
CORSIKA approach (section 2.2.2), where additionally the temperature and 
density dependence of the fluorescence yield is taken into account. 
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2.2.1 NKG approach 



In the usual treatment that was also used in a previous study of the shower 
image [5] , the fluorescence hght emitted by a shower is calculated from [2] : 

photon 



dldO: in 



sr m 



(2) 



where n^fi is a constant value of total fluorescence yield. The total number of 
particles Ne — J pN{X,r)2nrdr is given by the Gaisser-Hillas function [1] 

max Xo)/\ 

exp ((Xma. - X)/X) (3) 

where N^ax is the number of particles at shower maximum given by [14] 

= 0.7597 * 10^ (4) 

and p]^{X,r) is density of electrons in the shower given by the Nishimura- 
Kamata-Greisen (NKG) formula [15] 

o.(X r) - ^^(—Y-'(l + ' V-^-^^ r(4.5-.) 

rj, W ^'^TM^ 27rr(.)r(4.5-2.)' 

X is the atmospheric slant depth, Xq the depth of flrst interaction, X^ax the 
depth of shower maximum, A the hadronic interaction length in air (commonly 
fixed to a value of 70 g/cm^), s the shower age parameter (s = 1 at shower 
maximum) and tm the Moliere radius. 

The Moliere radius is a natural transverse scale set by multiple scattering, and 
it determines the lateral spread of the shower. Since the electron radiation 
length (the cascade unit) in air depends on temperature and pressure, the 
Mohere radius varies along the shower path. The distribution of particles in 
a shower at a given depth depends on the history of the changes of tm along 
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the shower path rather than on the local tm value at this depth. To take this 
into account, one uses the tm value calculated at 2 cascade units above the 
current depth [16]: 

J, jj^j ^ P[mb] -73.94 cos e -J 1./5.25588 

In the NKG approach we keep a constant value of fluorescence yield n^fl = 4.02 
photons per meter, as used by the HiRes group [2]. The spatial distribution 
of emitted light is therefore also given by the NKG formula, and the shape 
function follows from Eq. (5) as f^Kci''") = PN{X,r)/Ne{X). The fluorescence 
light fraction -Fnkg(^) using equation (1) can then be determined analytically 
by the normalized incomplete beta function, 

X 

i^NKG(r) = 4(a, h) = — i- / u'^-\l - uf-'du (7) 

B{a,h) J 

where x = 1/ {I+Tm/t), a = s,b = 4.5— 2s and B{a, b) is Euler's beta function. 
Using the series expansion of Ix{a,b), [17] the fluorescence light fraction can 
be given by 

1 1 / ^ B(s + l.n + l) ( 1 



For s — 1 formula (8) reduces to 

i^NKG(r) = 1 - (l + — ) (9) 



Inverting the above equation and taking into account the distance from the 
detector to the shower (i?) we can flnd the angular size q;nkg that corresponds 
to a certain fraction of the total fluorescence light signal: 

aNKG = 2arctan = 2 arctan(^((l - FNKG(r))-°-^ - 1)). (10) 
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2.2.2 CORSIKA approach based on energy deposit 



In contrast to the NKG approach, the fluorescence hght production in the 
CORSIKA approach is connected to the local energy release of the shower 
particles in the air; additionally, a dependence of the yield on the local atmo- 
spheric conditions is taken into account: 



A dE 



photon 



m 



(11) 



where eA(P, T) is the fluorescence efficiency; pair, P and T are density, pressure 
and temperature of air, respectively; A is the photon wavelength, c is speed of 
hght and h is the Planck constant. 

In the CORSIKA shower simulation program, the energy loss dE/dX of the 
shower particles is calculated in detail, taking into account also the contribu- 
tion of particles below the simulation energy threshold [7]. We extended the 
code to obtain a spatial distribution of the energy deposit. This offers the pos- 
sibility to construct a shower simulation chain which allows the comparison of 
quantities very close to the measured ones, e.g. photon flux or distribution of 
light received at the detector or even per pixel as a function of time. In par- 
ticular, shower-to-shower fluctuations generated by CORSIKA are preserved 
in this way. 

The adopted air shower simulation part of the simulation chain is illustrated 
in Figure 2. A two-dimensional energy deposit distribution around the shower 
axis is stored in histograms during the simulation process for different atmo- 
spheric depths. By interpolation between the difi^crcnt atmospheric levels, a 
complete description of the spatial energy deposit distribution of the shower, 
taking into account also the geometrical time delays, is achieved. More specifl- 
cally, the lateral energy deposit psi^i, r) is calculated for 20 horizontal layers 
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of AX = 1 g/cm^. Each observation level corresponds to a certain vertical 
atmospheric depth, the first one to Xi — 120 g/cm^ and the last one to 
X20 = 870 g/cm^. 

The simultaneous photons, which constitute an instantaneous image of the 
shower, originate from a range of shower development stages [5], from the 
surface S as shown in Figure 2. These simultaneous photons are defined as 
those which arrive at the FD during a short time window At. During this At 
(corresponding to a small change of the shower position in the sky by Ax = 
0.04° as chosen in the code) the shower front moves downward along the shower 
axis by a small distance -RA^. This means that the small element of surface S 
in polar coordinates corresponds to a small volume AV = r A0Ari?Ax, where 
A0 and Ar are steps in the azimuth angle and in the radial direction relative 
to the shower axis and R is the distance from the FD to the volume AV . 
The volume AV^ is located between two CORSIKA observation levels Xi and 
Xj+i. The distance between these two levels is divided into N sublevels, each 
of them labeled by n. Due to the small spacing between the chosen CORSIKA 
levels, the value of energy deposit within the volume AV at distance r can 
then be constructed sufficiently well by linear interpolation: 

(V \ - n)pE{Xi, r) + npEjXj+i, r) 
pE[Xn,r)^ — . (12) 

An additional linear interpolation in radial direction between bins of the COR- 
SIKA output was used to find the density Pb(X„, r) of the energy deposit. ^ 



^ The step used in radial direction is Ar = 1 m and the binning of the two- 
dimensional CORSIKA histograms of energy deposit is 1 m xl m at distances 
smaller than 20 m to shower axis, and lOmxlOmat larger distances. 
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Using the above interpolation, the number of photons Nj from the volume 
AV that are emitted towards the FD can be calculated as: 

where i runs over 16 wavelength bins, ej is the fluorescence yield for a 1.4 
MeV electron at pressure of 760 mm Hg and temperature of 14° C, dS is 
a projection of the surface rA0Ar into surface perpendicular to direction 
of the shower axis, {dE / dX)\iAMeV is the electron energy loss evaluated at 
1.4 MeV, A is the hght collecting area of the FD, Rp is the shower impact 
parameter with respect to the FD and gi{p,T) is a function describing the 
dependence of the fluorescence yield on the density p and temperature T of 
the air. Kakimoto et al. [11] provided an analytical formula for gi{p,T). For 
the 391.4 nm fluorescence line (13th bin in formula (13)) 

and for the rest of the fluorescence spectrum 

9i{p, T) = ^ (15) 

' 2.760Fi(l + p5iVT) ^ ' 

where p is in units of g/cvo? and T is in Kelvin. Fi, Ai, A21 Bi and B2 are 
constants and are 1.044xl0-^ 0.929 cm^g-^ 0.574 CTa^g'^, 1850 cm^g-^K'^/^ 
6500 cm^g~^K~^/^, respectively. The value of 2.760 photon/m is the total 
fluorescence yield outside the 391 nm band. 

In the CORSIKA simulations performed for this analysis, electromagnetic in- 
teractions are treated by an upgraded version [18] of the EGS4 [19] code. 
High-energy hadronic interactions are calculated by the QGSJET [20] inter- 
action model. To reduce computing time for the simulation of high-energy 
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events, a thinning algorithm [21] is selected within CORSIKA: Only a subset 
of the secondary particles that have energies below a specified fraction of the 
primary energy are tracked in detail. An appropriate weight is attached to 
each tracked particle to assure energy conservation. The artificial fluctuations 
introduced by thinning are sufficiently small, when a thinning level of 10~^ 
with optimum weight limitation [18,22,23] has been chosen. This weight lim- 
itation stops thinning in case of large particle weights and includes different 
weight limits for the electromagnetic component compared to the muonic and 
hadronic ones. 

3 Results 

Simulation runs were performed for proton and iron showers for different pri- 
mary energy Eq. The depth of first interaction Xq in the NKG approximation 
was chosen according to the average depth Xq from CORSIKA, see Table 1. 
Showers landing at variable core distance Rp = 2,3, ...,11 and 12 km were 
studied. The results shown in the following refer to the shower maximum, 
where also the fiuorescence emission is largest. 

3.1 Shower image in the NKG approach 

The shape function of particle density /nkg(?^) at shower maximum is shown 
in Figure 3 A for vertical and inclined = 45°) showers with £'o = 10 EeV and 
Eq = 100 EeV. It can be seen that these shape functions are almost identical. 
Some differences between vertical and inclined showers are seen only at dis- 
tances to shower axis larger than ~ 50 m. The differences are due to changes 
of the Moliere radius with altitude: a larger zenith angle of the shower implies 
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a higher position of the shower maximum and in consequence a larger value 
of the Moliere radius. Since the Moliere radius determines the lateral spread 
of particles in the shower, for inchned showers the shape function /nkgI'") 
becomes broader. A similar effect can be observed for showers with the same 
geometry, but different energies (showers with lower energy have a higher po- 
sition of the maximum and also a larger Mohere radius) but in these cases the 
differences are much smaller. 

In the NKG approach the size of the shower image q;nkg is connected to the 
width of the shape function of particle density /nkg('") a-nd can be calculated 
at shower maximum using Eq. (10) for fixed Mohere radius, fraction of fluores- 
cence light Fnkg(?^) and the detector-to-shower distance R. The appropriate 
-^nkg('") functions for showers presented in Figure 3A are shown in Figure SB. 
It is seen that 90% (67%) of fluorescence light emitted (i.e. of shower par- 
ticles) are found within distances about 160 m (58 m) around shower axis 
for vertical showers and about 190 m (70 m) for inclined showers. The cor- 
responding angular width of fluorescence light distributions at the detector, 
positioned for instance at R = 3.16 km, in these cases is about 5.7° (2°) for 
vertical showers and 7.0° (2.6°) for inchned showers. In Table 2 the sizes of 
shower image containing 90% and 67% of the signal according to formula (10) 
are listed. There is about 5% difference in the image spot size of showers with 
the same zenith angle but different energy, and about 19% between inclined 
and vertical showers. In Figure 3C the dependence of shower image versus R 
in the NKG approach is shown. The 90% spot size exceeds 1.5° for vertical 
(inclined) showers at distances smaller than 12 km (14.5 km). With typical 
FD pixel resolution of 1-1.5°, the shower image will cover more than one pixel 
at these distances. For a correct primary energy determination of these events, 
the ffuorescence light in the neighboring pixels has to be taken into account. 
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For example, at R — 3.16 km the fraction of light outside the circle corre- 
sponding to pixel field of view (1.5° in diameter) is about 40%, as marked by 
the vertical dashed line in Figure 3B, but only about 10% if the R increases 4 
times (increasing R leads to proportional decreasing of image size) . Neglecting 
this effect would result in a significant underestimation of the reconstructed 
primary energy, especially for nearby showers. The analysis of Figure 3 and 
Table 2 leads to the following conclusion: in case of the NKG approximation 
the size of the shower image is independent of the primary energy for show- 
ers at the same development stage and geometry. In other words, the same 
Moliere radius and shower age imply the same shape of /nkg('") function and 
in consequence lead to the same spot size of the shower image. 

In the above estimation of shower image we have neglected the influence of 
Rayleigh- and Mie-scattered and direct Cherenkov light distribution on the 
shower image size. To estimate this effect, relative differences between shower 
images obtained using these additional contributions to the fluorescence flux 
with respect to fluorescence only are shown in Figure 4. The additional contri- 
butions to the fluorescence light increase the image size on average by about 
7% (3%) within the image size containing 90% (67%) of light and these changes 
of shower image size slightly depend on R. These changes can be well under- 
stood if we take into account the Rayleigh scattering, which is the second 
dominant component in the total signal for the studied geometry. It is well 
known that Rayleigh scattering probabihty is proportional to (1 -|- cos^^), 
where ^ is the angle between the direction of photon emission and the direc- 
tion towards the FD. Since ^ increases for a vertical shower with increasing 
R, so the Rayleigh scattering probability and also contribution of Rayleigh- 
scattered hght to the shower image will be smaller. We note from Figure 4 that 
this contribution depends on the fraction of light considered: it is larger when 
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we study 90% of light than for 67%. This means that in the "center" of shower 
image fluorescence dominates, but it is less in the "tail". The shower image in 
the scattered light is therefore larger, although the "scattered" contribution 
is small for the considered geometries. In the following we concentrate on the 
main component, the unscattered fluorescence light. 

3.2 Comparison of shower image in the NKG and CORSIKA approaches 

In this section we study the differences between the calculated lateral distri- 
butions of energy deposit in the NKG and CORSIKA approaches and their 
influence on the shower image. We assume that fluorescence emmision dom- 
inates the received signal and that the distribution of light emitted by the 
shower is proportional to the distribution of energy deposit: /(r) ~ r). 
For this purpose, in Figure 5A we show the calculated lateral distributions 
of the energy deposit versus the distances to the shower axis at any point of 
surface S (see Figure 2). In case of the NKG approximation, the lateral density 
of energy deposit (dashed hne) is calculated using the following formula: 

PNKg(^, r) = {dE/dX)NmaxfNKG{r) (16) 

where (dE/dX) is the energy loss of an electron corresponding to a constant 
value of the average fluorescence yield n^^ — 4.02 photons per meter. 

In case of the CORSIKA approach, the energy deposit density (sohd line in 
Figure 5A) was obtained using the two-dimensional histogram of dE/dX. It is 
seen that the density of energy deposit obtained using CORSIKA histograms 
becomes larger than NKG at distances to shower axis smaller than 45 m. In 
the NKG approximation, it is assumed that all particles lose the same amount 
of energy and that the shape of the lateral distribution of energy deposit has 
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the same (NKG) functional form. Plots in Figures 5B and 5C show that these 
assumptions are not strictly valid. In Figure 5B we see that the particle density 
calculated from the NKG formula (dashed line) is different from the particle 
density from CORSIKA (solid line). The difference in the lateral distribution 
of energy deposit is mainly caused by this difference in the lateral particle 
distribution. A minor additional effect on the shape function is given by the 
average energy loss per particle. In Figure 5C the calculated relative difference 
z — 1 — {dE /dX) / {dE / dX) COR between average energy losses of electrons 
(dE/dX) in CORSIKA and NKG approach is shown. The average CORSIKA 
energy loss is always larger than energy loss in the NKG approach and the 
differences varies with distance from shower axis. This reflects a variation of 
the distribution of kinetic energy of particles around the shower axis, with 
more energetic particles being closer to the axis. Qualitatively, a narrower 
lateral particle distribution is expected for the CORSIKA proton events, as 
the electromagnetic component is permanently fed from high-energy hadrons 
coUimated around the axis. The NKG approximation, on the contrary, rather 
reflects a purely electromagnetic shower behavior. 

In Figure 6A, the normalized distribution of energy deposit from Figure 5 A 
(the shape function of energy deposit /^(r)) in the NKG and CORSIKA ap- 
proximations are shown. We see that for distance to shower axis smaller than 
25 m the CORSIKA shape function becomes considerably larger than the 
NKG one. Fitting CORSIKA data with a NKG-type function with fixed age 
s = 1 leads to an effective value of the Mohere radius = 58 m. This value 
is about 50% smaller than the original Moliere radius {tm = 104 m) in the 
NKG approach. This implies that the differences in the NKG and CORSIKA 
approaches will lead to different sizes of shower image. To estimate this dif- 
ference more precisely, first we calculate the fraction of energy deposit ^^(r) 
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based on /^(r) in CORSIKA and NKG approaches (see Figure 6B). Next we 
fit a two-parameter function 

FE{r)^l-{l + ^y\ (17) 

which is motivated by the functional form derived in Eq. (9), to the fraction 
of energy deposit. The fit leads to the following values of parameters a — 
54.24 ± 1.53 m and b = 1.928 ± 0.033. Using the above parameterization of 
Fe{'>^), we find the angular size qjcor corresponding to a given percentage of 
fiuorescence fight Fsir) in the CORSIKA approach: 

acoR = 2 arctan(|((l - FE(r))-'/' - 1)). (18) 

The size of shower image c^nkg in the NKG approach can be calculated using 
Eq. (10). In Figure 6C the shower image size ccnkg and cccor containing 90% of 
light are shown. We can see that the shower image in NKG approach is larger 
by about 23% than in CORSIKA. Finally, we calculate the relative difference 
k between the size of shower image in NKG and CORSIKA approach as a 
function of percentage of fiuorescence light according to the following formula: 



^ _ «NKG - QCOR ^ _ Q (1 - FEjr)) - 1 

"NKG ~ rjw (1 - -pNKG(r))-°-'^ - 1 ' 



(19) 



The variation of k is shown in Figure 6D. 

3.3 Shower image in the CORSIKA approach 

3.3.1 Dependence on primary energy 

The shape functions of CORSIKA lateral distributions for proton showers 
with primary energies £'o = 10 EeV (dashed line) and 100 EeV (solid fine) 
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are shown in Figure 7A. It is seen that the higher energy leads to a shghtly 
narrower shape function for distances above 10 m to shower axis. This imphes 
that the size of the shower image will decrease with increasing energy. Figures 7 
B, C and D confirm this result. The variation of the image size is rather small: 
below 7% in full FE{r) range. We note that the variation of the image size 
with energy is almost the same as that in the NKG approach (section 3.1). 



3.3.2 Dependence on zenith angle 

The integral of energy deposit Fsir) for proton showers with zenith angles 
^ = 0" and 6 = 45° at energy 10 EeV is shown in Figure 8A. 90% of the 
energy deposit is found within the distance of 125 m for ^ = 0" and 170 m 
for 9 — 45° around the shower axis. This means that the image spot size is 
about 4.52° and 6.15°, respectively (see also Table 3). A fit of a functional 
form as given in Eq. (17) to the fraction FE{r) of the energy deposit leads 
to = 137.2 ± 3.4 m and 645 = 2.86 ± 0.05 for the inchned shower and 
to oo — 54.24 ± 1.53 m and 60 — 1-928 ± 0.033 for the vertical shower. Using 
these parameters, we can find the angular size of the shower image according to 
formula (18) and the relative difference between showers with different zenith 
angles: 

0^45° - «o° gp (1 - F^(r))-V'"' - 1 

= " ' - ^(1-F.(r))-V^«-1- (2°) 

where a^o and a^^o are angular sizes of shower image for ^ = 0° and 6 — 
45°, respectively. The ratio /ccor versus fraction of hght FE{r) is shown in 
Figure 8B. It is interesting to compare these differences with those observed 
in the NKG approach. In the NKG approach, the size of the shower image 
depends on the Moliere radius (equation (10)), so for the same fraction of 
light FeIt) the relative differences for shower with different zenith angle is 
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given by 



^NKG — 



Q;NKG,45° — ttNKG,0° 
aNKG,45° 



r M,0° 



(21) 



where rM,ib° and rM,o° are the Mohere radii corresponding to the position 
of shower maximum for inchned and vertical shower, respectively. Using tm 
values from Table 2, we obtain /cnkg = 19%. We note that this value does 
not depend on the fraction of light -FB(r) (horizontal line in Figure 8B), in 
contrast to the difference /ccor, which strongly decreases with FE{r). 

3.3.3 Dependence on primary particle type 

Average lateral distributions of energy deposit in showers with different pri- 
mary particle and energy are presented in Figure 9. The lines represent three- 
parameter fits of NKG-type functions to data points; the parameters are shown 
in Table 4. The and s are only effective fitting parameters, not "real" 
Moliere radius and age parameter. The NKG function describes the COR- 
SIKA distribution of energy deposit very well close to shower axis, but with 
non-conventional and s. ^ It seems that such parameterization will be use- 
ful to calculate quickly the fluorescence signal using formula (13). Variation 
of the parameters (r^, s) with energy is not strong. For instance, in case of 
proton showers varies by about 2% between 10 and 100 EeV and the s 
parameter varies by about 9%. This means that at first approximation the 
shape of energy deposit density around the shower maximum seems to be 
almost independent of energy, although the amount of total energy deposit 
changes. On the other hand, when we compare s and for showers with the 
same energy but different primary, the differences are much larger. 
^ fitting with constant age parameters leads to worse jndf value, as shown in 



Table 4. 
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On the basis of Figure 9, one expects differences in the size of shower image for 
the same energy, but different primary. To study this effect more precisely, we 
show in Figure lOB the integral of energy deposit -FE(r) for iron and proton 
shape function at 10 EeV. It can be seen that 90% of energy deposit falls 
within 125 m from the shower axis in case of proton shower, and within 149 
m for iron shower. The image spot size is about 4.5° and 5.4° for proton 
and iron, respectively. A fit of a functional form as given by Eq. 17 to the 
fraction of energy deposit in iron showers leads to ape — 55.79 ± 1.83 m 
and bpe — 1-805 ± 0.038 and in proton showers ap — 54.24 ± 1.53 m and 
bp = 1.928 ± 0.033. Thus, the size of shower image for iron showers ape and 
proton one ap can be calculated using formula (18) with appropriate values 
of parameters; an example is shown in Figure IOC. The size of iron shower 
image is always larger by about 13% than proton one for all distances. In 
Figure lOD the relative difference q = {aFe — (^p)/(^Fe between iron and proton 
shower image size versus fraction of light is presented. It can be seen that 
differences in the image spot size between iron and proton increase when we 
take into account a larger fraction of the energy deposit. We note that the 
difference q was calculated assuming the same distance to the shower, but 
not the same altitude of the proton and iron shower maximum. It should 
therefore be checked if the observed difference between iron and proton image 
is only an atmospheric effect given by the different local value of tm in air- 
This atmospheric effect can be estimated using the Mohere radius for proton 
and iron showers at their maxima and can be calculated using the equation 
q = 1 — rM,Fe/rM,p- Since the Moliere radius for iron rM,Fe = 110 m and 
for proton rM,p = 104 m, the relative difference in the shower image due to 
the atmospheric effect is ^ ~ —6%. Thus, half of the difference between the 
primaries visible in q presented in Figure lOD is caused by this atmospheric 
effect. 
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In Figure 11 the influence of fluctuations in proton and iron shower shape 
function of energy deposit are presented. Fluctuations in proton shower profile 
lead to changes in the size of the shower image of about 1°. However, the image 
of a proton shower is always smaller than iron shower image. 

3.4 Detailed simulations of the shower image 

This section summarizes results presented until now with one modification: 
we show the shower image including all light components. 

Figure 12 shows the size of the shower image a containing 90% or 67% of 
light as a function of distance R from the FD to the shower, for showers with 
different core distance Rp. A comparison of the shower image derived using the 
two-dimensional CORSIKA histograms and that given by the NKG function 
is made for two different shower energies. It is evident that the image size 
in the shower maximum is independent of energy in the NKG approximation 
and that the NKG approximation leads to larger sizes of shower image than 
those derived from CORSIKA. Moreover, for a shower with higher energy the 
image size from CORSIKA is slightly smaller than the size at lower energy. 
These differences can be understood when we take into account the variation 
of the shape function in these cases, which were discussed earlier and shown 
in Figures 3A, 6A and 7A. For example. Figure 7A shows that the values of 
the shape function at 100 EeV are larger than those at 10 EeV at distances 
to the shower axis smaller than 10 m. Since we calculate the widths of these 
functions at distances corresponding to 90% or 67% of the total signal, wc 
expect that the width at the higher energy will be smaller. A similar effect is 
observed when one compares the shape functions in the NKG and CORSIKA 
approximations, (see Figure 6A). In this case one expects that the width of 
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the shape function in the NKG approximation will be larger than that derived 
from CORSIKA. In case of the NKG approximation, the changes of the shape 
functions with energy are negligible, as seen in Figure 3A; the observed small 
differences are only due to different distances to the shower. 

4 Conclusion 

Shower image simulations more accurate than available until now are pre- 
sented, which incorporate a more realistic distribution of fluorescence light 
emitted by the shower. The image simulations are based on distributions of 
energy deposited by the shower in air as derived from CORSIKA. A com- 
parison of the size of the shower image obtained using CORSIKA and that 
given by the NKG function was made for different energies and primary parti- 
cles. To a first approximation, the results of these two completely independent 
methods (analytical versus Monte Carlo) show quite reasonable agreement. 

The image spot size derived from CORSIKA is smaller by about 15% compared 
to the NKG approximation. This difference is mainly due to the differences in 
lateral particle distributions in the NKG and CORSIKA approximation. 

The energy deposit distribution from CORSIKA leads to a dependence of the 
size of shower image on the primary particle, so that studies of the shower 
image may be helpful for the primary particle identification. 
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Table 1 

Average values of depth of first interaction Xq, depth of shower maximum X^nax 
and altitude of shower max Hz (above sea level) for vertical showers obtained from 
CORSIKA. 



Eq Xq X^ax Hz 



(EeV) (g/cm^) (g/cm^) (km) 

p 10 44.4 757 2.572 

100 42.1 805 2.034 

Fe 10 10.6 696 3.241 

100 8.7 746 2.695 
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Table 2 

Size of shower image a and distance r around the shower axis containing 90% and 
67% of fluorescence light in the NKG approximation for vertical proton showers of 
different Eq and zenith angle 6 landing at a distance of 3 km from the eye and 
observed from R=3.16 km. Additionally, the local Moliere radius tm is shown at 
shower maxima. 

^0 090% ^^90% «67% ^^67% I'm 

(EeV) (deg) (dcg) (m) (deg) (m) (m) 



10 5.69 157 2.10 58 104 

10 45 7.00 194 2.59 71 128 

100 5.42 150 2.00 55 99 

100 45 6.68 184 2.47 68 122 



Table 3 

Size of shower image a and distance r containig 90% and 67% of fluorescence light 
for vertical proton showers of different zenith angle 6 at energy 10 EeV. 

«90% ''90% "67% ^67% 

(deg) (deg) (m) (deg) (m) 



4.52 125 1.53 42 
45 6.15 170 2.26 65 
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Table 4 

Fitting parameters of NKG-type functions (r„t, s and Nmax at shower maximum) 
to shape functions obtained using CORSIKA lateral distribution of energy deposit. 
^m(s = 1) is the value obtained using fixed values of s parameter at shower maxi- 



mum. 



Eq Vm s Nmax / ndf r^(s=l) x^/ndf 

(EeV) (m) (10^° particles) (m) 



10 98.1 ±0.2 0.844 ±0.001 1.572 ±0.001 3.95 58 4.1 

100 96.7 ±0.2 0.765 ±0.001 16.492 ± 0.004 2.00 46 8.3 



Fe 10 46.5 ±0.9 1.181 ± 0.009 1.532 ±0.001 1.26 68 1.6 

100 46.6 ±0.8 1.201 ±0.008 15.349 ± 0.008 1.22 65 1.7 
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Fig. 1. Energy dependence of nitrogen fluorescence between 300 and 400 nm in dry 
air at the pressure 760 mm Hg. The scale of fluorescence yield is adjusted so that 
the 1.4 MeV point lies on the dE/dX curve (taken from Ref. [llj)- 




Fig. 2. Geometry of an EAS as seen by the fluorescence detector. Photons which 
arrive simultaneously to the FD originate from surface S. See text for more details. 
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Fig. 3. (A) Shape of particle density distribution /nkg('') versus distance to shower 
axis in the NKG approximation. Showers with different energies Eq and zenith an- 
gles 6 are shown; (B) Integral -Fnkg('') of the shape function /nkg('") from Fig- 
ure 3 A. Horizontal dashed lines correspond to 90% and 67%. The vertical dashed 
line indicates a 1.5° pixel detector field of view. Upper scale is the shower image 
size corresponding to distance to shower R=3.16 km. (C) Size of shower image 
containing 90% of fluorescence light versus the detector-to- shower distance (R). 
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Fig. 4. Relative difference between size of shower image calculated using total light 
and only fluorescence versus the detector-to-shower (R) distance for a vertical 
shower with energy 10 EeV. 
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Fig. 5. (A) Lateral distributions of energy deposit density in the CORSIKA and 
NKG approximations, calculated for an average vertical proton shower with energy 
10 EeV. (B) Particle density from CORSIKA and derived using NKG function, (C) 
Relative difference z between average energy loss obtained from the CORSIKA and 
NKG approaches. 
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Fig. 6. (A) Shape of particles density in the CORSIKA and NKG approx- 

imations. (B) Integral FE{r) of shape functions fEir) from Figure 6A. (C) Size 
of shower image containing 90% of fluorescence light versus the detector-to- shower 
distance R. (D) Relative difference between shower image size obtained in the COR- 
SIKA and NKC approaches, see text for more details. Vertical showers at energy 10 
EeV are presented. 
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Fig. 7. (A) Shape of energy deposit density fE(r) for vertical proton showers with 
energies 100 EcV and 10 EcV derived from CORSIKA. (B) Integral of energy de- 
posits versus distance to shower axis. (C) Size of shower image containing 90% of 
fluorescence light versus detector-to-shower distance R. (D) Relative difference in 
the shower image size between proton showers with energies 10 EeV and 100 EeV. 
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Fig. 8. (A) Integral of energy deposits versus distance to shower axis for proton 
showers with different inclination, derived from CORSIKA. (B) Relative difference 
in the shower size image between these showers. 
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Fig. 9. Comparison of the the average lateral distributions of energy deposit density 
calculated for an average iron and proton showers with different energy. 
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Fig. 10. (A) Shape of energy deposit density fE{r) for vertical 10 EeV proton and 

iron showers. (B) Integral of energy deposit versus distance to shower axis for proton 
shower (solid line) and iron (dashed line). (C) Size of shower image containing 90% 
of fluorescence light versus detector to shower distance R. (D) Relative difference 
in the shower image between iron and proton shower. 
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Fig. 11. (A) Lateral distribution of energy deposit calculated for 15 single proton 
showers (solid lines) and average of 5 iron showers (dashed lines) at energy Eq = 10 
EeV. (B) Integral of the energy deposits versus distance to shower axis for proton 
showers (solid lines) and iron ones (dashed line). 
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Fig. 12. (A) Size of the shower image at shower maximum containing 90% and 
67% of light versus the detector to shower distance R, using the CORSIKA and 
NKG distributions of energy deposit. The dashed line corresponds to shower image 
obtained with constant value of fluorescence yield n^fi = 4.02 photons/m. (B) Size 
of the shower image containing 90% and 67% of light versus R using the CORSIKA 
distributions of energy deposit for iron (showers at different energies). 
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